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We compute the dimensionally regularised four-loop vacuum energy density of the SU(-/V C ) 
gauge + adjoint Higgs theory, in the disordered phase. "Scalarisation" , or reduction to a 
small set of master integrals of the type appearing in scalar field theories, is carried out in d 
dimensions, employing general partial integration identities through an algorithm developed 
by Laporta, while the remaining scalar integrals are evaluated in d = 3 — 2e dimensions, 
by expanding in e < 1 and evaluating a number of coefficients. The results have implica- 
tions for the thermodynamics of finite temperature QCD, allowing to determine perturbative 
contributions of orders 0(g G ln(l/g)), 0(g e ) to the pressure, while the general methods are 
applicable also to studies of critical phenomena in QED-like statistical physics systems. 
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1. Introduction 



The theory we study in this paper is the Euclidean SU(iV c ) gauge + adjoint Higgs theory, 
defined in continuum dimensional regularisation by the action 

Se = Jd d x£E, (1.1) 

C E = l^Fii +Tr[D k ,A } 2 + m 2 TtA 2 + \(TrA 2 ) 2 , (1.2) 

where k,l = 1, ...,d, D k = 8 k - igA k , A k = A a k T a , A = A a Q T a , F kl = (i/g)[D k ,Di], and T a 
are the Hermitean generators of SU(iV c ), normalised as Tr T a T b = 5 ab /2. Summation over 
repeated indices is understood. We could have taken the scalar potential also in the form 
Ai(Tr Aq) 2 + A2Tr^Q, but the two quartic terms are independent only for N c > 4 and thus, 
to avoid further proliferation of formulae, we will set A2 = here, denoting A = Ai. For the 
moment we keep d general, but later on we write d = 3 — 2e, and expand in e <C 1. 

The observable we would like to compute for the theory in Eq. (|1.2|) is its partition function, 
or "vacuum energy density", 

f(m 2 ,g 2 ,X) = - lim i In [vA k VA exp(-S E ) . (1.3) 

Here V is the d-dimensional volume. The phase diagram of the system described by Se 
has a "disordered", or symmetric phase and, depending on N c , various kinds of symmetry 
broken phases ^ [2] ■ Our aim is to determine the perturbative expansion for / up to 4- loop 
order in the symmetric phase, expanding around Aq = Af, = 0; the 3-loop result is known 
already [31 II]- The result will depend on N c through dA = N 2 — 1,Ca = Nc- 

The main motivation for the exercise described comes from finite temperature QCD. Indeed, 
the simplest physical observable there, the free energy density or minus the pressure, has been 
computed perturbatively up to resummed 3-loop level jSJ |S] , but the expansion converges 
very slowly, requiring probably temperatures T TeV to make any sense at all [3J El El 
|7]. Moreover, at the 4-loop level the expansion breaks down completely [HI El- Multiloop 
computations are not useless, though: these infrared problems can be isolated into the three- 
dimensional (3d) effective field theory in Eq. (|1.2|) [TO] , and studied non-perturbatively there 
with simple lattice simulations 4 . However, to convert the results from lattice regularisation 
to 3d continuum regularisation, and from the 3d continuum theory to the original 4d physical 
theory, still necessitates a number of perturbative multiloop "matching" computations. 

The way our computation enters this setup has been described in Combining our re- 
sults with those of another paper ^2] allows one to determine, as explained in jllj . all the log- 
arithmic ultraviolet and infrared divergences entering the 4-loop free energy of QCD. This not 
only fixes the last perturbatively computable contribution to the free energy of hot QCD |llj . 
of order 0(g e ln(l/g)T 4 ), but is also a step towards renormalising the non-perturbative con- 
tributions, as determined with lattice methods [31ll3|. Some other applications of our results 
are discussed in Sec. |H1 
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Figure 1: The skeleton diagrams contributing in Eq. after subtraction of those which 

obviously vanish because of colour contractions or specific properties of dimensional regulari- 
sation. Solid lines represent the adjoint scalar Aq, wavy lines the gauge boson A{, and dotted 
lines the ghosts. The complete sets of skeleton diagrams have been enumerated and written 



2. Outline of the general procedure 

The first step of the perturbative computation is the generation of the Feynman diagrams. 
At 4-loop level, this is no longer a completely trivial task. In order to make the procedure 
tractable, we employ an algorithm whereby the graphs are generated in two sets: two-particle- 
irreducible "skeleton" graphs, as well as various types of "ring" diagrams, containing all 
possible self-energy insertions. The resulting sets, with the relevant symmetry factors, were 
provided explicitly in |14j . 

It actually turns out that some of the generic graphs shown in ^3] do not contribute in 
the present computation. There are two reasons for this. First, once the Feynman rules 
for the interactions of gauge bosons and adjoint scalars are taken into account, some of the 
graphs vanish at the point of colour contractions. This concerns particularly the "non-planar" 
topologies [IK]. Second, all vacuum graphs which do not contain at least one massive (adjoint 
scalar) line, vanish in strict dimensional regularisation. In some cases such a vanishing may 
be due to an unphysical cancellation between ultraviolet and infrared divergences, as we will 
recall in Sec.0 but for the moment we accept the vanishing literally. The remaining skeleton 
graphs are then as shown in Fig. ^ For the ring diagrams, which by far outnumber the 
skeleton graphs, we find it simpler to treat the full sets as shown in ^3], letting the two 
types of cancellations mentioned above come out automatically in the actual computation. 
For completeness, the ring diagrams are reproduced in Fig. El 

The Feynman rules for the vertices and propagators appearing are the standard ones. We 



down in refi whose overall sign conventions we also follow. 
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Figure 2: T/ie ring diagrams contributing in Eq. kl.tyl ' L lJfi. The notation 
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employ covariant gauge fixing, with a general gauge fixing parameter, denoted here by 

£ = Chore = 1 ^standard > (^-l) 

where Standard is the gauge fixing parameter of the standard covariant gauges. Therefore, 
Feynman gauge corresponds here to £ = 0, Landau gauge to £ = 1. We keep everywhere £ 
completely general, however, and verify explicitly that it cancels in all the results. 

The graphs having been identified and the Feynman rules specified, we program them 
in the symbolic manipulation package FORM |17j . for further treatment. 

After the colour contractions, the next step is to "scalarise" the remaining integrals. That 
is, we want to remove all scalar products from the numerators of the momentum integrations, 
such that only integrations of the type appearing in scalar field theories remain. This problem 
can be solved by using general partial integration identities |18| . The full power of the 
identities can be conveniently made use of through an algorithm developed by Laporta |19j . 
We discuss some aspects of our implementation of this algorithm, together with the results 
obtained, in Sec. El 

After the reduction to scalar integrals, we are faced with their evaluation. At this point one 
has to specify the dimension d of the spacetime, in order to make further progress. We write 
d = 3 — 2e, expand in e <C 1, and evaluate the various scalar integrals appearing to a certain 
(integral-dependent) depth in this expansion, such that a specified order is achieved for the 
overall result. For the new 4-loop contributions, the overall order for which we have either 
analytic or numerical expressions is 0(1). The scalar integrals needed for this are discussed 
in Sec. 0] 

There is one remaining step to be taken before we have the final result: the renormalisation 
of the parameters m 2 ,g 2 , A in Eq. ()1.2[) . In other words, the results presented up to this point 
were in terms of the bare parameters, and we now want to re-expand them in terms of the 
renormalised parameters. This step is also specific to the dimension, and turns out to be 
particularly simple for d = 3 — 2e, since only the mass parameter gets renormalised. The 
conversion of the bare parameters to the renormalised ones is discussed in Sec. El and the 
final form of the results is then shown in Sec. El 

Having completed the straightforward computation, we discuss the conceptual issue of 
infrared divergences in Sec. We mention in this context also some checks of our results, 
based on largely independent computations. We end with a list of some applications in Sec. [HJ 

3. Scalarisation in d dimensions 

After inserting the Feynman rules and carrying out the colour contractions, there remains, at 
4-loop level, a 4<i-dimensional momentum integration to be carried out. The different types of 
integrations emerging can be illustrated in graphical notation in the standard way. Without 
specifying the fairly complicated numerators, involving all possible kinds of scalar products 
of the integration momenta, the graphs are of the general types shown in Fig. El 



4 



(H> ; @> (M> 

(S> (3D' O*0 ; OQZ> 
O*0> 0*0; (O) 2 *0; (O) 4 

Figure 3: The 15 general types of 4-loop integrations remaining, in terms of momentum flow 
(momentum conservation is assumed at the vertices), after taking into account that colour 
contractions remove the non-planar topologies. Any line could contain a propagator to some 
power n > 1, and there is also an unspecified collection of scalar products of the integration 
momenta in the numerator. 



There are a few simple tricks available in order to try and simplify the scalar products 
appearing in the numerators 16 . For instance, one can find relabelings of the integration 
variables such that the denominators appearing in the graph remain the same, while the 
scalar products in the numerators may get simplified, after symmetrising between such rela- 
belings. Some scalar products in the numerators can also be completed into sums of squares, 
such that they cancel against the denominators. Furthermore, we can make use of various 
special properties of dimensional regularisation: any closed massless 1-loop tadpole integral 
vanishes; and any 1-loop massive bubble diagram with at most one external momentum is 
easily scalarised explicitly, in the sense of removing the loop momentum from all the scalar 
products appearing in the numerators. However, while such simple tricks are sufficient at, 
say, 2-loop level, this is no longer the case at 4-loop level. 

To scalarise the 4-loop integrations, we have to make full use of the identities provided by 
general partial integrations ^Sl- To systematically employ all such identities, we implement 
the algorithm presented by Laporta ^S] using the "tables" routines of FORM |17j . This leads 
to a complete solution of our problem. The main technical details of our implementation were 
discussed in [TB] . 

After the scalarisation, the master integrals remaining are those shown in Fig. |1J This 
basis is, of course, not unique. As an example, one could have chosen a different basis for the 
3-loop master integrals, employing identities following from partial integrations |20| . 
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Figure 4: T/ie 1-loop, 3-loop and J^-loop "master" topologies remaining after "scalarisation" . 
There are no numerators left in these graphs. A solid line is a massive propagator, l/(p 2 + 
m 2 ), and a dotted line a massless one, 1/p 2 , where p is the Euclidean momentum flowing 
through the line. Note that no independent 2-loop representative appears. 



and correspondingly for the higher loop integrals. Therefore, the 3-loop master integrals we 
are using, appearing on the right-hand-sides of Eqs. ()3.1j) . (|3.2j) . could be exchanged in favour 
of the 3-loop integrals on the left-hand-sides of Eqs. ()3.1|) . (|3.2I) . 

To display the full result after scalarisation, we introduce the shorthand notations 



f(m 2 ,g 2 ,X) 



A = 



A(d A + 2) 
9 2 C A 



We then obtain the following expressions: 
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a 2 = -14081760 + 11237380d + 64451424d 2 - 140115669d 3 + 129957772d 4 

-69456108d 5 + 23323366d 6 - 5020699d 7 + 674926d 8 - 51720d 9 + 1728d 
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It is worth stressing that Eqs. (|3.5|) - (|3.14j) were obtained with an arbitrary £, which just 
exactly cancelled once all the graphs were summed together, for a general d, and before 
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inserting any properties of the master integrals. This is a consequence of the fact that the 
master integrals constitute really a linearly independent basis for the present problem. 



4. Integrals in d = 3 — 2e dimensions 



A set of master scalar integrals having been identified, the next step is to compute them. As 
already mentioned, we do this by writing d = 3 — 2e, expanding in e 1, and evaluating a 
number of coefficients in the series. 

In order to display the results, we first choose a convenient integration measure. To this 
end, we introduce an MS scale parameter p,, by writing each integration as 

id m r 2e (eiy f d d p 
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where /i = /i(e 7 /47r) 1//2 , and the expression in square brackets has integer dimensionality. 
This square bracket part of an £-loop integration is then written as 
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where counts the number of propagators, or lines, in the graphical representation of the 
function g. From now on we assume that the loop integrations are computed with the 
dimensionless measure in the curly brackets in Eq. (|4.2[l . while the constants in front of the 
curly brackets, together with the explicit powers of m as they appear in Eqs. (|3.5|)H|3.8[) . are 
to be provided in trivial prefactors (cf. Eq. (|4.14j) below). 

With such conventions, the loop integrals remaining are functions of e only, and read: 
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As we will see, the terms shown explicitly are needed for determining the 1/e-poles in the 
4- loop expression for /, the constants 7„ are needed for determining the finite 4- loop contri- 
bution to /, and the higher order terms only contribute at the level 0(e). Analytic results 
for 71, ...,79, as well as a numerical determination of 710, are presented in |Appendix A 



It is now convenient to combine the conventions in Eqs. (|3.4j) . ()4.1j) . ()4.2|) and write 
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Substituting Eqs. (|4.3f) — (|4. 13f) into Eqs. (|3.5j) ~ H3.8j) and expanding in e, the results then read, 
up to 0(e) corrections: 
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It is interesting to note that while single diagrams contributing to p3 do have 1/e-poles (cf. 
Appendix B ), they sum to zero in the term without A, but not in the terms proportional to 



A, A 2 . This structure is related to counterterm contributions from lower orders, as discussed in 
the next section. Similarly, single diagrams contributing to p4 have both 1/e 2 and 1/e-poles, 
but the former ones sum to zero in the term without any A's. 

Of course, single diagrams contain also ^-dependence. In our computation £ cancelled at 
the stage of Eqs. ()3, 5j) - (j3.14|) . but one could alternatively express single diagrams in terms 
of the same basis of master integrals, this time with ^-dependent coefficients, and let the £'s 
sum to zero only in the end. For completeness, we again illustrate the general structure of 
such expressions at the 3-loop level, in Appendix B 



5. Counterterm contributions 

The computation so far has been in terms of the bare parameters of the Lagrangian in 
Eq. (|1.2|) . As a final step the result is, however, to be converted into an expansion in terms 
of the renormalised parameters. 

The conversion is particularly simple in low dimensions such as close to d = 3, since 
then the theory in Eq. (J1.2|) is super-renormalisable. In fact, the only parameter requiring 
renormalisation is the mass parameter m 2 . We write it as 

m 2 = m 2 are = m 2 (fi) + 5m 2 , (5-1) 

5m 2 = 2{d A + 2)j—^^-(-g 2 \C A + \ 2 ). (5.2) 

This exact counterterm |21| I22j guarantees that all n-point Green's functions computed with 
the theory are ultraviolet finite. Note that as far as dimensional reasons and single diagrams 
are concerned, there could also be divergences of the form g 4 /e, but they sum to zero in the 
counterterm appearing in Eq. (|5.2I) . 

Inserting now Eqs. (|5,1|) . (|5.2|) into the 1-loop and 2-loop expressions for f(m 2 ,g 2 ,X), we 
get contributions of the same order as the 3-loop and 4- loop vacuum graphs, respectively, 
from 5m 2 ■ d m 2 f (m 2 , g 2 , A). We need to use here Eqs. (|3.5|) . (|3.6() . since 0(e)-terms, not 
shown in Eqs. ([4.15ft . (|4.16j) . contribute as well, being multiplied by the 1/e in 5m 2 . Explicitly, 
the terms to be added to Eqs. (|4.17[) . (|4.18j) . once the prefactors in Eq. (|4.14j) are expressed 
in terms of the renormalised parameter m(/2) rather than m, are 



I » V'f 1 , !V < , 1 



* - U + 5A-* + 4T+5*J- (5 ' 3) 

* - (24))" 4, B)(?^ (1+x,+ 5 (4 ^ 2)+2S )(-^^ 2 )- (M) 

The 3-loop 1/e-contributions in Eq. ()5.3|) cancel against the 1/e-terms in Eq. (|4.17jl . Indeed, 
genuine vacuum divergences can only appear in p2,P4, since such divergences must be analytic 
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in the parameters m ,g ,X appearing in the Lagrangian, while p% comes with a coefficient 
~ (m 2 (/x)) 1//2 (cf. Eq. I|4.14|) ). Another point to note is that l/e 2 -terms appear in 5p^ only 
with coefficients A, A 2 , just as in Eq. (|4,18|) . although there is no complete cancellation. 



6. The final result 

We can now collect together the full result for f(m 2 ,g 2 ,X), in terms of the renormalised 
parameters of the theory. For dimensional reasons, its structure is, 



-2e r 



47T 



+ 



+ 



+ 



-4e r 



(4vr) 2 



(47T) 



(47T) 



/l,0 m 3 (j2) + 
/2,0 5 2 + / 2 ,iA]m 2 (/Z) + 

/3,0 9 4 + /3,i5 2 A + /3, 2 A 2 ]m(/Z) + 
/ 4 ,o / + /4,i /A + / 4 , 2 5 2 A 2 + / 4l3 A E 



+ 



(6.1) 



where fa = fi i(e, /J,/m(fi)) are dimensionless numbers, with £ indicating the loop order, and 
i the number of A's appearing: 



/i,o 

/2,0 
hi 
/3,0 
/3,1 

/3,2 
/4,0 



d A {d A + 2)[- +0(e) 
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7T _ 



d ACi[ — - — ln2 + — +C(e 



24 6 



6 



d A C A (d A + 2)(-\n-^—--+0{e, 
V lm{pL) 4 



cU(cU + 2)(lm 



A* 



2m(/z) 2 
+ d A (d A + 2) 2 (-i +C?(e)) , 



ln2 + 0(e) + 



(6.2) 
(6.3) 
(6.4) 
(6.5) 
(6.6) 

(6.7) 



d A C\ 



43 491 



32 6144 
85291 1259 



-7T 



•--81*-^ 



+ 



In 2 



5653 



1 



19 



1536" +4^" 3 72 + y 73 + 



768 32 

157 13. . 4 19 1 . 

+777^4 + ^t(75 + 76 ) + - 77 + t^78 + t 79 - 7io + O(e) 



/4,1 



192 

<UC 2 



64 



64 



-7T 



1 



48 ' 6 
p 



(6. 



8 In 



12 



- H In 2 vr 2 + - (74 + 75) + 0(e) 

2 2 192 4 w ' ; v ; 



+ d A C\{d A + 2) 



+ 



16e 2 



+ In 2 — ^— + - In ■ 

2m(/2) 2 2m(/2) 



24 12 12 v 7 



(6.9) 



/4,2 



^^(^ + 2) 



32 



7T 



16e 2 



64e 



In 2 



183 51 , 

H n 2 

4 2 



23 
48 



1 



2m(/i) 
1 



/36-vr 2 
V 8 



In 



2m(/u) 



vr 2 + 273 + - 74 + - 7 8 + C(e 



+ 



/ 



1,3 



+ d A C A ( ( i A + 2) 2 (^ In 

= d A (d A + 2)(d A + 

+ d A (d A + 2) 2 (-^ In 

+ d A (d A + 2) 



1 

2 

,3(1 

V24 



1 

— In 

2 2m(/2) 

V192 
[1 



(6.10) 



+ 81n o , 

^ +Iln2 + 0(e; 



2m(jtt) 
+ 0(e 

In particular, following the notation of ref. jllj and writing 



g 74 + O(e) 



+ 



+ 



/4,0 = — G? A C A 



nxj, ; + 81n 2^) 



we read from Eq. Q6.8|) that 
43 491 

«M = 
/?M = 



32 6144 
85291 1259 



-vr 2 ss 0.555017 



(6.11) 



(6.12) 



(6.13) 



768 
157 

74 

192 ' 



32 
13 



In 2 + 



5653 



-7T 



1 5 

I 7i + g 72 

1 



19 



311 
256 
-1.391512 . 



1536 

. 4 19 
64 (75 + 76) - 9 77 " 48 78 

19, 2 77 2 
— In 2 ln 2 2 + vr 2 

32 6 9216 1536 



-73 



1. 



79 + 7io 



43 



491 2 , 1793 

vr 2 In 2 + 

512 



C(3)+7io 



(6.14) 



In Eq. (|6.14|) we used values for 71, ...,710 from Appendix IA.5l and Appendix IA. 61 The coef- 
ficient «m gives a contribution of order 0(<? 6 ln(l/<?)T 4 ) and @m a perturbative contribution 
of order 0(g 6 T 4 ) to the pressure of hot QCD jUJ. 

The expression in Eq. (|6.1I) . with the coefficients in Eqs. H6.2|) - ()6.11|) . contains a number of 
1/e 2 and 1/e-poles. Once our computation is embedded into some physical setting, such as 
in a vacuum counterterm is automatically generated (denoted by pe(T) in ^J), which 
eventually cancels all the UV-poles, such that physical observables remain finite for e — > 0. 
The nature of the poles in Eqs. (|6.2|) - Q6.11|) is analysed in detail in the next section. 
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7. Infrared insensitivity of the results 

The result shown in the previous section contains a number of 1/e 2 and 1/e-divergences. Since 
dimensional regularisation regulates at the same time both ultraviolet (UV) and infrared (IR) 
divergences, we may ask of what type are those obtained? The purpose of this section is to 
show that the divergences are of purely UV origin, and the result is thus IR insensitive, if 
interpreted properly. There are two ways of showing this, firstly an effective theory approach 
in which one understands that all the IR divergences are contained in the SU(iV c ) pure Yang- 
Mills theory obtained by integrating out the Ao-field, secondly a pragmatic one in which one 
shields away the IR divergences by giving the gluon and ghost fields a mass. 

Conceptually the best way to analyse the IR sensitivity is to dress the problem in an 
effective theory language. In the present context, such an analysis was carried out in |1U| . 
The idea is that since the field Aq has a mass scale, it can be integrated out. The integration 
out is an ultraviolet procedure, thus by construction not sensitive to IR physics. The effective 
low-energy theory that emerges is a 3d pure gauge theory. Its partition function, on the other 
hand, does contain IR divergences, starting at 4-loop level jHUHj- 

Therefore, we expect that all results up to 3-loop level should be IR insensitive. At 4-loop 
level there is a part of the result, that is the diagrams which can be constructed fully inside 
the pure SU(iV c ) theory, which can be both IR and UV divergent. Since in dimensional 
regularisation, however, these graphs are set to zero, the non-zero result we have obtained 
should again be insensitive to any mass scales in the gluon and ghost propagators. 

Apart from the issue mentioned, there is also another possible source of IR problems, 
namely that of overlapping divergences. Indeed, while IR divergences appear for vacuum 
graphs at 4-loop level only, they appear for self-energy graphs already at the 2-loop level 
(see, e.g., [TH]). However, 2-loop self-energy insertions do appear also as subgraphs in the 
4-loop "ring diagrams", making the divergence structure of such 4-loop graphs "doubly" 
problematic. We return to this issue presently, but first finish the discussion of IR divergences 
at lower than 4-loop level. 

To be very explicit, let us introduce a fictitious mass parameter trq for all massless lines 
(gluons and ghosts), hence giving the function f(m 2 ,g 2 ,X) a further functional dependence 
on the mass ratio x = ms/m. Let us denote by AH ("Adjoint Higgs") graphs with at least 
one ^o-line, and by YM ( "Yang- Mills" ) graphs with none at all. The general structure of the 
bare f(m 2 ,g 2 ,X) can then be expressed as (cf. Eq. (|4.14[) ) 

/<»v.^£(^) f (£) 2 'V)<-^- 

where /^ H , /^ YM are dimensionless functions. While the treatment above corresponds to 
setting x = first and then computing the expansion in e, we now keep a non-zero x through 
the entire calculation, being interested in the limit of small x only in the end: 

MS: limhm/A H (x,e,e,A) , (7.2) 



f £ Ati (x,e,Z,\) + x 



// M ( e ,0 , (7.1) 
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IR-regulator : lim lim f^ n (x, e, £, A) . (7.3) 

These two limits do not in general commute for single diagrams, but should commute for the 
sum. Possible power IR divergences in single diagrams would show up as poles in x, while 
logarithmic ones correspond to lncc. 

The main technical differences in the IR regularised procedure with respect to the MS 
computation are a more complicated scalarisation, in the absence of low-level routines specific 
to the presence of massless lines, such as the so-called "triangle rule", and an enlarged set 
of master integrals. Furthermore, some additional diagrams contribute, which were set to 
zero from the outset in the MS calculation, due to the absence of any mass scale (in some 
subdiagram). 

As a roundup, it turns out that, starting at the 3-loop level, individual diagrams do indeed 
contain logarithmic as well as powerlike IR divergences, which then cancel in the sum, proving 
a posteriori the validity of the dimensionally regularised MS calculation. For completeness, 



we illustrate this issue in Appendix B 



We now return to the 4-loop level. According to the discussion above, the full set of graphs 
can be divided into four sub-classes, having potentially different IR properties: pure Yang- 
Mills graphs (YM) and those with at least one ^o-line (AH), with both sets further divided 
into skeletons (Fig. ^) and ring diagrams (Fig. [2J). The properties of the pure Yang-Mills 
diagrams are discussed in ^5], and we only state here that they contain both logarithmic UV 
as well as IR divergences, which however exactly cancel in strict dimensional regularisation 
(but not in regularisations which only regulate the UV, such as lattice regularisation). Here 
we then just discuss the skeletons and rings containing at least one massive ^o-line. F° r 
simplicity, we discuss explicitly only terms without a quartic coupling A. 

We have computed the 1 /e-divergence in the sum of such AH-skeletons with in total three 
different mass spectra: 

1. As described above, whereby the Ao-lines carry the mass parameter m 2 , while the gluon 
and ghost lines are massless. 

2. By giving an equal mass to all the fields: Aq, gluons, and ghosts. The computation 
proceeds in complete analogy with the one described in |12j . 



3. By setting all masses to zero, picking some line in the 4-loop vacuum graph, integrat- 
ing the massless 3-loop 2-point function connected to that line in d dimensions 5 , and 
regulating the remaining single integral by shielding the IR with a mass and regulating 
the UV via dimensional regularization. 

All three methods give the same result for the 1/e-pole in AH-skeletons, confirming its ex- 
pected IR finiteness. Expressed as a contribution to p4 in Eq. (|4,18|) . the divergence appearing 

J This problem has been solved a long time ago via integration by parts; for a discussion as well as an 
algorithmic implementation, see |2M| . 
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in the result reads 



3072e 



1 




For the AH-rings, on the other hand, the third method does not work. This is due to the 
overlapping divergences mentioned above: a 2-loop 2-point function of gluons alone leads to 
logarithmic UV and IR divergences, and trying to carry out the final integration by some 
recipe, gives generically an outcome ~ 1/e 2 , but with a coefficient dependent on what the 
recipe precisely was. To cancel the l/e 2 -pole, not to mention to get the correct coefficient 
for the remaining 1/e-pole, is a very delicate problem, which can only be guaranteed to have 
been solved by employing a fully systematic procedure. Our non-Abelian case is therefore 
qualitatively different from a pure scalar theory, where the problem of overlapping divergences 
does not emerge |24| . For a discussion of the cancellation of the analogue of the l/e 2 -pole in 
cutoff regularisation in the pure SU(iV c ) theory, see [13] . 

On the contrary, the AH-rings can be systematically computed with the 1st and 2nd types of 
mass spectra. Both procedures give the same result, confirming its IR insensitivity. Summing 
together with Eq. (|7.4j) . we recover the ^-independent 1/e-pole on the first row in Eq. (|4.18j) . 

In summary, we have verified explicitly that the only possible IR divergence appearing in 
our computation is that of the pure SXJ(N C ) gauge theory, contained in the YM-graphs. It is 
addressed further in ref. |12j . 

8. Discussion and conclusions 

The main point of this paper has been the discussion of formal analytic techniques for, and 
actual results from, the evaluation of the 4-loop partition function of the 3d S\J(N C ) + adjoint 
Higgs theory using dimensional regularisation. The final result is shown in Eqs. ()6.1j) - (|6,llj) . 
We have also demonstrated that if interpreted as a matching coefficient — that is, if the 
pure Yang-Mills graphs, without any adjoint scalar lines, are dropped, as is automatically 
the case in strict dimensional regularisation — then the result is IR finite. Therefore, all 
IR divergences are contained in the pure Yang-Mills theory. We would now like to end by 
recalling that such techniques and results have also practical applications. 

Perhaps the most important application is that our results provide two specific new per- 
turbative contributions to the free energy of hot QCD, of orders g® ln(l / g)T 4 , g 6 T 4 [TT] . 
Similarly, they provide also new perturbative contributions to quark number susceptibili- 
ties |25| . Once the parameters of the 3d theory are expressed in terms of the parameters 
of the physical finite temperature QCD via dimensional reduction, and once other contribu- 
tions of the same parametric magnitudes are added, this allows for instance to re-estimate the 
convergence properties of QCD perturbation theory at high temperatures |11 ( I25[ 126]. Our 
present computation also contributes to the MS scheme renormalisation of the simplest 3d 
gauge-invariant local condensates, obtained by partial derivatives of the action with respect 
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to various parameters |U, an d thus in principle helps in non-perturbative studies of the 
pressure of high-temperature QCD [IJ^j. It may also allow for refined analytic estimates 
such as Pade resummations |27j for the observable in Eq. 1)1.3(1 . 

Let us mention that there has recently been significant interest in somewhat more phe- 
nomenological approaches to QCD perturbation theory at high temperatures (for reviews 
see, e.g., |2S|)- As far as we can tell our results are of no immediate use in such settings. 

Another application is that our general procedure is relevant for studies of critical phe- 
nomena in some statistical physics systems. In this context one may either study directly 
the three-dimensional physical system, or carry out computations first in d = 4 — e dimen- 
sions, expand in e, and then take the limit e — > 1. For instance, some properties of the 
Ginzburg-Landau theory of superconductivity have been addressed in the former setup up 
to 2-loop level (see, e.g., [2Hl _ [II2])j an d in the latter setup, in the disordered phase, up to 
3- loop level [HS1- The integrals arising in the disordered phase are "QED-like" just as in our 
study, so that scalarisation and the sets of master integrals are essentially the same as the 
present ones |16) . Moreover, in the case d = 4 — e, the master integrals can be evaluated to 
a high accuracy utilising the techniques introduced in |34j . while for d = 3 — 2e most master 
integrals have been evaluated in this paper. Our methods could therefore help in reaching 
the 4-loop order. 
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Appendix A. Master integrals 

We discuss in this Appendix the determination of the scalar master integrals of Eqs. 1)4.3(1 - 
1)4.13(1 . They depend on one mass-scale m only and are thus "QED-like" in the generalised 
sense that the number of massive lines at each vertex is even. Since the dependence on 
m is trivial and has been absorbed into the coefficients, see Eq. 1)4.2(1 . m = 1 in most of 
what follows. One obtains particularly simple expansions in 3 — 2e dimensions by using the 
integration measure / = Jd 3 ~ 2t p, in accordance with Eq. (|4.2j) . 

We first discuss briefly the various general techniques we have employed for the evaluation 
of these integrals. The list of techniques includes: partial integration relations between var- 
ious scalar integrals, in analogy with those derived at 3-loop level in j^Uj ( Appendix IA.1J) ; 
graphs with only two massive lines, which can often be evaluated exactly (Appendix IA.2() ; 
graphs with two or three vertices, which can be evaluated to a sufficient depth in e using con- 
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figuration space methods ( Appendix I A. 3|) : and some remaining graphs, which were evaluated 
in momentum space ( App endix I A . 4|) . We combine the results from the various techniques in 
Appendix lA.5| showing the actual expansions for the master integrals to the depths specified 
in Eqs. (|4.3|) - (|4.11|) . There is one finite integral remaining which we have not been able 
to evaluate analytically, corresponding to Eq. (|4.12j) : its numerical value is determined in 
Appendix IA. 61 



A.l. Partial integration identities 



Implementing systematically all identities following from partial integrations, as discussed 
in Sec. |HJ allows not only to express all integrals in terms of a few scalar ones, which do not 
contain any non-trivial numerators, but produces also a set of relations between the scalar 
integrals. As a simple example, we may recall that the identity 







±f — 

Jp,q 9pk 



Pk 



(p 2 + m 2 )(q 2 + m 2 )(p — q)' 



leads to the relation 



P,q (P 2 + m 2 )(q 2 + m 2 )(p — q) 2 d — 3 J p (p 2 + m 2 ) 2 J q (q 2 + m 2 ) 



(A.l) 



(A.2) 



Taking furthermore into account that in dimensional regularisation the two integrals on the 
right-hand-side of Eq. ()A.2|> are related, we obtain 



O m 2 (Q 



(d-2) 



2(d-3)j- (A - 3) 

Examples of similar relations at 3- loop level were shown in Eqs. (|3.1|) . (|3.2|) . and a complete 
3- loop analysis can be found in [2U] (see also jH]). 

At 4-loop level, there are obviously many more identities than at 3-loop level. Rather than 
showing a complete list we give here, as an example, one of the relations: 



d 



2d 



4(d - 3) 



(A.4) 



2(d - 3) 

It turns out that this relation is convenient for the determination of the 4-loop integral on 
the right-hand-side. 



A.2. Integrals known exactly 

A few of the integrals appearing can be evaluated exactly in d dimensions. This holds 
particularly for cases where only two massive propagators appear. As an example, we show 
how this can be done in configuration space. The massive propagator can be written as 

G i:m; = / - — r-5— ^ — k tt = r, — — Ki ArriiX) , (A. 5) 

V ' l! J (27r) 3 ~ 2 ' p 2 + m 2 (2vr)|- eV z/ 2 - eV ; v ' 
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where K is a modified Bessel function, and x denotes, depending on the context, either a 
d-dimensional vector or its modulus. On the other hand, the massless part of the graph 
converts in configuration space to 



( 2vr )3-2e 



1 



r(|) 



2 i/ 7T2 ~ e X 



£™3— 2e 



We can then employ the identity 



} A-2 



dx x x K 2 (x) 



r(i + A)- 



1 + A + 2/i\ 2 /l + A\/l + A - 2ijl 



(A.6) 



(A.7) 



With this result, the following expressions are easily derived (using the integration measure 
inside the curly brackets in Eq. (|4.2j) ): 
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Obviously we also know ( f is again according to the curly brackets in Eq. (|4.2|l ) 
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(A.8) 
(A.9) 
(A.10) 
(A.ll) 
(A.12) 
(A.13) 

(A.14) 
(A.15) 



A. 3. Configuration space evaluations 

Even when configuration space does not allow for an exact evaluation of the integral, like 
in Appendix IA.21 it may allow for the most straightforward way of obtaining a number of 
coefficients in an expansion of the result in e. This is the case particularly if there are only 
two vertices in the graph. 

At £-loop level, the graphs in this class are of the form 
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where G(x;mi) is from Eq. (|A.5f) . The idea (see, e.g., |24| ^ is to split the integration into 
two parts: J °°dx(...) = J r dx(...) + f r °° dx(...). The first part is performed in d = 3 — 2e 
dimensions but by using the asymptotic small- X form of G(x;m,i), 



G(x; rrii 



r(j-c) i 

4vrl- e ^" 2e 



U J r(|- e) + 



mix 



+ .. 



(A.17) 



while the latter part, which is finite, is performed by expanding first in e and then carrying 
out the remaining integrals is d = 3 dimensions. For instance, 



G(x;rrii) 
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+ C(e 2 



. (A.18) 



When the two parts are summed together and the limit r — > is taken, the dependence on r 
cancels, and we obtain the desired result. 

In the evaluation of such integrals, dilogarithms will in general appear. Their properties 
have been summarised, e.g., in [33]. For completeness, let us recall here that one can shift 
the argument of 



Li 2 (x) 



n>0 U 



(A.19) 



from the intervals [-00...-I], [-1...0], [1/2. ..1] to the interval [0...1/2] via 
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respectively. As follows from here, the dilogarithms satisfy, for x > 0, 
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(A.20) 
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Special values include 

Li 2 (-1) = -— , Li 2 (0) = 0, Li 2 (-)= _ . ... , , 

12 ' v ; ' V2 / 12 2 ' 6 

Using these identities, and denoting M = mi + m 2 + 7713, we obtain for the 2- loop case 
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For the 3-loop case, now denoting M = mi + 7712 + 7773 + 777,4, we obtain 
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In particular, if all masses are equal, 

=™ -8 + 41n2-4e(l3 + ^vr 2 -81n2 + ln 2 2) + 0(e 
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The case of two massless and two massive lines can be checked against Eq. I)A,9|) . 

The 4-loop case has only been worked out to order 0(1), rather than 0(e). Denoting now 

M = mi + 7772 + m% + 777-4 + TU5, 
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(A.28) 



The case of three massless and two massive lines can be checked against Eq. (|A.13|) . The 
next order, 0(e), could also be worked out and is indeed needed for 7§ in Eq. (|4.1U[) . but we 
choose to use another way to determine it, based on Eq. (|A.4[) . 

When there are more than two vertices in the graph, the configuration space technique 
gets rapidly more complicated, due to the difficult structure of the angular integrals. There 
is one graph we are interested in, however, whose divergent and, most incredibly, also the 
constant part ^Hi can still be obtained analytically: 



2 3\ /4 5 6 



47T 



(?)■ 



d : 



3-2e 



X 



3 6 

d^y G(x - y; mi) J[ G(x; m*) J] G(y; mj) . (A.29) 

i=2 j=4 



Employing the angular integral [37] 
K x (\x-y\) (2*)™ 



1 



l x — y\ x 



(xy)' 



9{x - y)K x {x)I x {y) + 9{y - x)K x (y)I x (x) 



(A.30) 



where A = \ — e and on the right-hand-side x = |z|, y = \y\, one is left with two inde- 
pendent radial integrations which can be handled as above [30]) by splitting the integra- 
tions as J °°dx(...) = J r dx(...) + J* r °° da; (...)• Denoting M i2 3 = mi + m 2 +7773, M23456 = 
7772 + 1^3 + m 4 + m 5 + Tt6> the outcome is 

1 2 \ 8e 1 



2 3\ /4 5 6 



MfflJ 32 



VM23456 M123 M123 



, (A.31) 
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where 



In particular, 



+ 



13 H vr 2 - 41n 2 x + 

12 



+ 2Li 2 (l -y) + 2Li 2 (y + z) + 2Li 2 (-z) +8 



1 — x 
x(l + z) 



Li 2 (l -x) + 



1 + 



1 - X 
x{l + z) 



TT- 



Li 2 (— xz) + lnx ln(l + xz) 

6 



1 

32 



i+-+4(l3-81n 2 2-^ 2 )+0(e 
e z e V 12 



(A.32) 



(A.33) 



A. 4. Momentum space evaluations 

When the graph has more than two vertices, the configuration space method is in general no 
longer practical. Some of these graphs are, however, rather easily evaluated in momentum 
space. This is the case particularly for the "triangle" topology, shown in Eq. (|A.38|) below. 
The triangle graph consists of three consecutive 1-loop self-energy insertions, 



d 3-2e g 



i r(i +6) 



B(p,m 1 ,m 2 ,e) 



(A.34) 



where B(p,mi,m2,c) is a one-dimensional integral over a Feynman parameter. It has the 
properties 



B(0, mi,m 2 ,e) = 0, 

B(p, 0, 0, e) = ]im^B(p,rni,rri2,c) 

P 



B(p, mi, ni2, 0) = 2arctan 



r(l-2e) : 



mi + fn 2 



(A.35) 
(A.36) 
(A.37) 



The triangle graph is then just a one-dimensional integration over the modulus of p. Carrying 
out one partial integration and expanding in e, one obtains |3S] 



8e TT Z 



where 



mi + m 2 / 32 



64 f 00 d 
X(x,y) = -? / dp hip — 
vr- 4 Jo dp 



- +2 + 41n2-x 
e 



/m 3 + m 4 m 5 + m 6 
\mi + m 2 ' mi + m 2 



d 
dp 



arctan(p) arctan I — I arctan I — 



V 



y 



OA ££j 

X (1,1) = ^C(3), x(l,0) = 4C(3), X(0,0)=0. 
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7T- 



(A.38) 

(A.39) 
(A.40) 
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A. 5. Summary of expansions for master integrals 



Given the results of the previous sections, we can collect together the expressions for the 
constants 71,. ..,79 defined in Eqs. (JO jH|4~TTjl . From Eq. ()A.15j) . 

71 = -8-y +|C(3). (A.41) 
Combining Eq. (|3~T|) with Eq. ()A.9|) . 

72 = g " \ C(3) • (A ' 42) 
Combining Eq. (|3~2~)) with Eq. (|A.27j) . 

73 = -^7r 2 -ln 2 2. (A.43) 

b 

From Eqs. (|A.38|) . (|A.40|) (or, for 7 6 , from Eq. (jA.lljO . 

^(2 + 41n2-^C(3)) , (A.44) 

^(2 + 41n2-^C(3)) , (A.45) 
2 

^(2 + 4 In 2), (A.46) 



74 

75 
76 



^(2 + 121n2-^C(3)) • (A.47) 



Combining Eqs. (|A.33[) . (fO|) . 

78 : 

Finally, from Eq. (|A.13j) . 



53 

175 - 96 In 2 + 16 ln 2 2 + — vr 2 . (A.48) 



33 55 

79 = 270 + __7r 2 -__C(3). (A.49) 



A. 6. Numerical evaluation of 7 10 

It can easily be verified that the integrals in Eqs. (|4.12l) . (|4. 13f) are both infrared and ultra- 
violet finite. They can therefore be evaluated directly in d = 3 dimensions. For the present 
application we only need 710, defined by Eq. (|4.12|) . 

There is no obvious partial integration relation whereby 710 could be reduced to a simpler 
integral. Due to the fact that there are four vertices, it is also not easily treated in configura- 
tion space. The most straightforward approach seems then to be to combine the self-energy 
of Eqs. (|A.34|) . (|A.37j) with the 2-loop self-energy 

(A.50) 
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e m co> gz> m gz> ©; © 



Figure 5: A possible choice for 1-loop, 2-loop and 3-loop "master" topologies, in the case that 
gluons and ghosts are treated as particles with a mass niQ. There are no numerators left in 
these graphs. A solid line is a propagator of the form l/(p 2 + m 2 ), and a dashed line of the 
form l/(p 2 +m 2 G ), where p is the Euclidean momentum flowing through the line. A line with 
a blob on it indicates a squared propagator, l/(p 2 +m 2 ) 2 . 



for which a one-dimensional integral representation has been given in |38| . This leads to a 
simple two-dimensional integral representation: 



7io 



2 f°° p 

— dp p arctan - • Tl2(p) , 

vr Jo 2 



(A.51) 



where (3S] 

n 2 (p) = 



i r l 



dx 



P S Jx-(p) y/p 2 X 2 - (1 - X 2 ) 2 I 1 - X 



P 



[1 + x 2 ) arctan 



P 



-2(1 — x + x 2 ) arctan 



l + x 2 V ' 



+ x In 



1 + 



P 



(1 + xf 



(A.52) 



and X-(p) = (l+p 2 /4) 2 —p/2. We may note that in Eq. ()A.52|) it is numerically advantageous 
to change the integration variable from x to y = y/x — X-(p). The final result reads 



710 ~ 0.171007009753(1) , 
where the number in parentheses indicates the uncertainty in the last digit. 



(A.53) 



Appendix B. Three-loop results with and without an IR cutoff 

As discussed in Seed starting at the 3-loop level single graphs are considerably more infrared 
sensitive than the total sum: the limits in Eqs. (|7.2j) . H7.3|) commute only for the latter. Let 
us recall that there x = mc/m, where niQ is a fictitious mass given to the gluons and ghosts: 

sab / c \ sab 

(A a k (p)AK-p)) = "if— -r (Ski -JW-aX-r) . (c a (p)c b (p)) = , (B.l) 

p z + tjiq V P + ttiq / p + m,Q 

where A^,c a ,c b are the gluon, ghost, and anti-ghost fields, respectively. We illustrate the 
general structures appearing here with a few specific examples. 
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In the presence of the two mass scales m, mo, the set of master integrals is more complicated 
than when gluons and ghosts are massless. The master integrals that can appear in principle 
are shown in Fig. up to 3-loop level. 

It turns out that the skeleton diagrams are better behaved in the IR than the ring diagrams: 
power and logarithmic IR divergences appear only in single rings, but they cancel in their 
sum. A rather typical example, with both an UV pole 1/e and an IR divergence ln(mG/m), 
is given by the gluon ring with a scalar and ghost bubble attached to it. Carrying out 
scalarisation to the master integrals shown in Fig. [5J denoting x = mc/m, and normalising 
as in Eq. (|3.7|) . we obtain 

(d-2) 



+ 

+ 
+ 

+ 



m 2 
m 2 



8(d- l)x 2 
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432 (d - 1) x 2 
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1134 - 459d + 54d 2 + 108£ - 



CD 



1296 (d - 1) 

90d£ - 18d 3 C - 84£ 2 + 76df + 5d 2 f + 2d 3 £ 2 + d 4 £ 2 ) + 
(32 - 12d - 6x 2 + 3dx 2 ) ' 



+ m 



32(d-l)x 2 
-1 + x 2 ) 



+ 



+ 



+ m 



2(d-l)x 2 _ 
(44 -8d- Ux 2 + 5dx 2 ) 



32 (d -i) j- < R2 » 

According to Eq. (j7.3j) . the first step is now to expand in e 1. The integrals emerging are 
all known [HJ • Changing the normalisation to be according to Eq. (|4.17f) , we obtain 
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1 3r 
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The second step is then to expand in x C 1: 
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+ 0(e) 



(B.3) 
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+ 0(x,e). (B.4) 
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We observe that there is a gauge-parameter dependent UV-divergence in the form of 1/e, and 
a gauge-parameter dependent logarithmic IR divergence in the form of lnx. 
Proceeding according to Eq. (|7.2|) . on the other hand, leads to 



i.d-2) 



+ m 



16(d-3)(2d-7)(3d 
(d-3) 



+ 



_4(2d-7)(3d-8) 
in terms of the master integrals in Fig. 0] Expanding in e <C 1 



3lEq. El 



3^n + ^ 



(B.5) 



(B.6) 



Clearly Eqs. (IB~4l . (lB~6l) do not agree 6 . Summing all the graphs together, however, both 
procedures lead to the gauge-parameter independent and UV and IR finite p3 on the first 
row in Eq. (|4.17jl : in other words, £, 1/e and lnx all cancel. 

Some other rings lead also to 1/x-divergences. Let us show, as an example, 



while 



5p 



3lEq. El 



1 24 
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12£ + 5£ 2 



64 



A +1 ln2 + 0(x,e) 



dp 



x=0 



3lEq. El 



— +- In2 + C(e) 
24 3 w 



(B.7) 



(B.8) 



Again, the 1/x-divergences of the type in Eq. (|B.7|) cancel when gluon rings with all possible 
1-loop scalar insertions are summed together. 

As a comparison of Eqs. ()B.4|) and (|B.6|) . or Eqs. ()B.7|) and (|B.8|I shows, the computa- 
tion carried out with an IR cutoff leads in general to a more pronounced gauge-parameter 
dependence for single graphs than the computation carried out according to Eq. (|7.2|) . just 
because the introduction of a mass according to Eq. ([B.ljl breaks gauge invariance. The 
results of Eqs. (jB.6|) . (|B.8|) are anomalously simple, however: in general there is certainly 
gauge-parameter dependence left over in single graphs also with the procedure of Eq. (j7.2J) . 
For example, 



Sp3 1 



x=0 



1 M3 if l & l 4 °- 2 ^ + 5 ^ 

Eq.E3-^( 13 -^ + ^ ) + ^ 



+ 0(e) , (B.9) 



and £ cancels only in the sum. 



6 The first two terms in Eq. fBlj can be written as (2 - 2£ + ( 2 )[\n(?,x /2) - l/(4e)]/8 + l/(32e), showing 
that the result of Eq. (IB.6II arises after a cancellation of IR and UV divergences in dimensional regularisation. 
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